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I. INTRODUCTION 


Astrophysics and neutrino physics entered a new era 
after the discovery of high-energy cosmic neutrino events 
observed by the IceCube Collaboration [l|,[3- The follow¬ 
up reports with additional statistics uncover the 

spectrum of cosmic neutrino in the energy range between 
0(10) TeV and 0(1) PeV. The reported spectrum shows 
some remarkable features, for example, (i) the neutrino 
flux diminishes steeply as the energy increases, and the 
best-fit spectral index is = 2.5 and (ii) there 
is a rap in the energy range between 400 TeV and 1 
PeV [ 1 , The high spectral index is consequent on 
the sudden end of the observed event spectrum at the 
high-energy edge {E,^ ~ 3 PeV) and high event rate at 
the low-energy bin < 100 TeV M- In contrast, in 
Refs. (Tol - [l^ . it is also pointed out that neutrino spec¬ 
trum from standard hadronuclear process {pp inelastic 
scattering) with the spectral index higher than s^, > 2.2 
causes a serious conflict with the gamma-ray observa¬ 
tion at Fermi-LAT, if the spectrum does not have cutoff 
at the low energy. The presence of the gap may also 
urge reconsideration of the assumption of simple power- 
law spectrum, which typically results from pp reaction 
in the so-called cosmic-ray reservoir such as the galaxy 
cluster. Although these features in the spectrum have 
not been conclusive in statistics, plenty of attempts have 
been made to reproduce them from the aspects of both 

• [l3) par- 


astrophysics 
tide physics 


(for a review, see Ref. 
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In Refs. [13-113: origin of the gap in the observed 

spectrum was asked to the attenuation of cosmic neu¬ 
trino, which is caused by the scattering with cosmic 
neutrino background (Cz/B) through a new interaction 
between neutrinos, the so-called neutrino secret inter¬ 
action [sO - Isljl PI In such a scenario, the narrow width 
of the gap can be explained by the resonant behaviour 
of the scattering. In the previous study [ 11 ], we intro¬ 
duced a new gauged leptonic force to explain the gap and 
pointed out that the leptonic force could simultaneously 
explain the disagreement between theory and experiment 
in muon anomalous magnetic moment. We improve in 
this paper our numerical method and calculate diffuse 
neutrino flux, taking account of the distribution of the 
source of cosmic neutrino with respect to the redshift. 
Moreover, we search through the model parameter space 
to find a set of parameters that can reproduce not only 
the gap but also the sharp edge at the upper end of the 
cosmic neutrino spectrum. The existence of the edge is 
expected to improve the fit of the spectrum with a lower 
value of spectral index to the observation. Here we also 
discuss constraints on the model, some of which were not 
considered in our previous study, such as the neutrino- 
electron scattering process, invisible decay of a light par¬ 
ticle at colliders, big bang nucleosynthesis (BBN), and 
supernova cooling. 

The paper is organized as follows: In the second sec¬ 
tion, we describe our model and illustrate parameter 


^ The resonant neutrino-neutrino scattering mediated by the stan¬ 
dard model (SM) Z boson has been investigated in the context 
of the Z-burst scenario to address ultrahigh-energy cosmic neu¬ 
trinos [sSlblH . For early works on the relation between the secret 
interaction and new mediation fields, see, e.g., Refs. [6H4fi9l| . The 
effect on the cosmic neutrino spectrum from the scatte ring be¬ 
tween neutrino and dark matter is discussed in Refs. [Ja. l47l. ITOII . 
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regions relevant to the cosmic neutrino spectrum and 
muon anomalous magnetic moment. We also discuss 
constraints from laboratory experiments and cosmologi¬ 
cal and astrophysical observations. Differential equations 
for diffuse neutrino flux with the leptonic force are given 
in Sec. uni The spectra calculated with different model 
parameters and redshift distribution of cosmic neutrino 
sources are compared with the observation in Sec. lIVI Fi¬ 
nally, we discuss the relation between the characteristic 
features of the cosmic neutrino spectrum and neutrino 
mass spectrum, and also the distribution of sources of 
cosmic neutrino. 



II. MODEL AND CONSTRAINTS 

We extend the SM of particle physics with a massive 
vector boson Z' that mediates a new leptonic force, 

•^int = gz'QaP [LaYLp + t-RaY^Rp] Z'p, (1) 

where and lua are a lepton doublet and a right- 
handed charged lepton singlet with flavour a = {e, p,, r}, 
respectively. We choose the flavour structure of the in¬ 
teraction as Qa^ = diag(0,1, —1), which corresponds to 
the 17(1) gauge interaction associated with muon number 
minus tau number [L^ — Lr) 0122. In this paper, we 
do not discuss the details of the model, such as a mecha¬ 
nism of the symmetry breaking. 0 Instead, we handle the 
two parameters, the coupling gz> and the mass Mz' of 
the gauge boson, as parameters that describe the model. 

The interaction with neutrinos in Eq. © is expected 
to produce the gap and the edge in the cosmic neutrino 
spectrum through the resonant scattering with Ci^B. In 
the resonant scattering process J^Cosmicf'Ci/B —>■ medi¬ 
ated by Z', only cosmic neutrinos having the energy cor¬ 
responding to the resonance energy ifres are selectively 
scattered off by Gi/B on the way from its source to the 
IceCube poMdll . 1^ 1^. which results in the gap around 
Eres- Here, E^es is given as 

^—- ( 2 ) 

2m^(l -I- z) 

where stands for a mass of the target CvB and 2 
is the redshift parameter at which the scattering occurs. 
In Eq. ([2]), the C^B is assumed to be at rest. With an 
assumption of rrii, = 0(0.1) eV, the scale of Mz' can 
be estimated as Mz' = 0(1 — 10) MeV for E^es — 1 
PeV. Meanwhile, in order to scatter a sufficient amount 
of cosmic neutrino during the travel of 0(1) Gpc, the size 
of the cross section is required to be larger than ^ 10“^° 


^ Although we do not fully describe the model Lagrangian, we 
assume that the Yukawa sector also respects the symmetry, and it 
is broken so as not to shift the mass eigenbasis of charged leptons, 
i.e., the lepton flavour is not violated. The gauge symmetry and 
its phenomenology have been discussed in Refs. I73ll82ll . 


FIG. 1. Relevant parameter region and the constraints. The 
red band represents a parameter region favoured by muon 
anomalous magnetic moment within 2a. The hatched region 
is excluded by the lepton trident search at the CCFR experi¬ 
ment. The region excluded by the measurement of ve —>■ ve 
at Borexino is filled with a gray colour. The vertical dashed 
line stands for the lower bound on Mz' from BBN. Two sym¬ 
bols + and X indicate (Mz', Qz') = (H MeV, 5 x 10“^) and 
(9 MeV, 4 X 10“^), respectively, which are used in Sec. IV. 
See the text for details. 


cm^ at the resonance [12 |4l| . In the — Lt model, the 
cross section near the resonance is estimated ai0 


^res 


M|, V s 


( 3 ) 


by using a d-function approximation, where s ~ 2m^Ei, 
is the square of the center-of-mass energy in the limit 
of small GvB momentum. The requirement to the cross 
section turns out to be gz' ^ 0(10“^). Putting it all 
together, the model parameter region that is relevant to 
the cosmic neutrino spectrum at the energy range around 
1 PeV can be deduced as 


gz' > 0(10"^) and Mz' = 0(1 ~ 10) MeV. (4) 

The interaction with charged leptons (namely muons) 
in Eq. o is responsible for an extra contribution 
to muon anomalous magnetic moment. The measure¬ 
ment [11] deviates from the SM predictions by 

around Str. The extra contributions from various types 
of new physics to fill this discrepancy have been inten¬ 
sively studied [11] , and also in the context of the — Lr 
model [12 • The Z' contribution with the combinations 
of gz' and Mz' indicated with the red band in Fig. [T] suc¬ 
cessfully reproduces the observed value of muon anoma¬ 
lous magnetic moment within 2cr errors. In Fig. [I] we 
also show some experimental constraints (see below) on 
the model, which pin down the model parameter region 
onto 


gz' 0(10-^) at Mz' = 0(10) MeV. (5) 


^ In the numerical calculation, the scattering cross section in the 
neutrino mass eigenbasis is used, cf. Eqs. 0 and ®. 
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It is worth noting that this parameter region has some 
overlap with the region Eq. (|31) required to reproduce the 
gap and the edge in the cosmic neutrino spectrum. 

In the following, we summarize experimental con¬ 
straints on the — Lt model. 


o Neutrino trident production process; In the mass 
region {Mz' < 100 MeV) we are working with, the tight¬ 
est constraint comes from searches for the neutrino tri¬ 
dent process: the production process of a pair with 

a muon neutrino, which results from the scattering of a 
muon neutrino from a target nucleus. The process was 
observed at several neutrino beam experiments, e.g., the 
CHARM-II and the CCFR 91|. Since the reported 
results are in good agreement with the SM prediction, 
which is mediated by the Z and the W bosons, an addi¬ 
tional contribution to the process mediated by Z', whose 
amplitude is proportional to g ^,, is strongly constrained. 
Indeed, the constraint on Mz' and gz' was recently eval¬ 
uated in Ref. [92|. We adopt the 95% C.L. bound of 
the CCFR experiment, which is shown in Fig. [T] as the 
hatched region. 


o Neutrino-electron scattering: Although the elec¬ 
tron is not charged under the — Lr symmetry, Z' 
can interact with an electron through a kinetic mixing e 
between Z' and a photon, which is induced by loop dia¬ 
grams. The total contribution to the kinetic mixing e in 
the — Lr model is hnite, and it is estimated as 


|cioop I 


8 egz' 
3 (47r)2 


In — = 7.2 • 10"® 

rrifj, 


( 


9Z' 

5 • 10-4 


). 


( 6 ) 


where e is the elementary electric charge. This leads 
to an extra contribution to the elastic ve —> ve scat¬ 
tering signal in the solar neutrino measurement at the 
Borexino experiment. The amplitude is proportional to 
eegz' — Mzr), where is the momentum transfer. 
The constraints from Borexino are discussed in Refs, 

1^ in the context of various different scenarios of new 
leptonic force. Here, we interpret the bounds given in 
Ref. [ 9 ^ with those in the — Lr model while taking 
account of the fraction of mass eigenstates of solar neu¬ 
trino, which is given in Ref. [^. In Fig. [1] the excluded 
region is filled with a gray colour. The measurement of v- 
e elastic scattering at LSND places a similar bound 
to the coupling gz' at Mzr ~ 10 MeV (and a weaker 
bound at Mz' — 1 MeV) 

o Beam dump experiments: Once the kinetic mixing 
between a photon and Z' appears, Z' can be produced 
on shell at beam dump experiments, whose production 
rate is O(e^e^). However, the Z' immediately decays into 
the neutrino-antineutrino channel in the — Lr model. 
Therefore, the constraints from the search for e^e” pairs 
at the detectors located downstream of the beam dump, 
such as E137 are irrelevant for our scenario. 

o Invisible decay of a light particle: Since Z' in 
this framework decays dominantly to one of the 
smoking-gun signals in collider experiments is a light par¬ 


ticle decaying to an invisible mode. The kinetic mix¬ 
ing leads the on-shell production of Z'. The pr ocess , 
e+e” —>• 'jZ' —>■ 7 -Finvisible, is searched at BABAR [l0r)j |. 
which sets the bound to the coupling between the elec¬ 
tron and Z' at ~ 10“^ |l0l| . Our reference choices of 
parameters satisfy this condition. 

o BBN: A constraint on Mz' is derived from BBN. If Z' 
is as light as the temperature at the era of BBN, its exis¬ 
tence increases the number of relativistic degrees of free¬ 
dom, Nes, and the success of the standard BBN might 
be sp oiled, which leads the lower bound Mz' > I MeV 
[I02j | . This condition is always satisfied on the parameter 
region of our interest. Nevertheless, Z' with a mass of 
0(10) MeV may indirectly contribute to Nes through a 
raise in the temperature of and iZr [^ . In Fig. [U we 
display the lower bound on Mz' from the indirect con¬ 
tribution with ATVeff < 0.7 as the vertical dashed line, 
which is taken from Ref. 

o Other constraints and future improvement: 

There exist many other experimental observations con¬ 
straining the Lf^ — Lr model. Most of them, however, are 
either irrelevant in the mass region we are focusing on 
(i.e. Mz' < 100 MeV) or weaker than the bounds dis¬ 
cussed above. Here, we mention some of them; searches 
for the SM Z boson decay to four leptons at the LHC 
can be used to constrain the — Lr model [zi [lii- 
However, it becomes insensitive in the mass region below 
Mz' — 10 GeV due to the invariant mass cut for the 
same-flavour leptons. Furthermore, Z' alters the decay 
rates of Z, W, and mesons, but constraint s from them 
are weaker than the CCFR bound [Tol fTo^. The bound 
from a precise measurement of cosmic microwave back¬ 
ground anisotropy is much weaker than those discussed 
above, cf. Fig. 1 in [4l| and discussion in Refs. fl05l. [l0^ 


It should be noted that a severe constraint on gz' pos¬ 
sibly arises from supernova neutrino observations |46l| . 
The existence of Z' with an MeV-scale mass interrupts 
free streaming of and Vr in the core, which could make 
the diffusion time much longer than the estimation from 
the standard supernova cooling process. To avoid possi¬ 
ble problems, an introduction of an additional invisible 
particle, e.g. the QCD axion, might be necessary. 


Apart from astrophysical/cosmological observations 
and the constraints through the kinetic mixing term, a 
direct test of a muonic force mediated by a bos on with a 
sub-GeV mass is planned at GERN SPS [l07j| . which is 
expected to improve the sensitivity to the coupling gz' by 
orders of magnitude and fully cover the parameter region 
referred with Eqs. o and ®. We expect that the tri¬ 
dent events have been recorded at near detectors in mod¬ 
ern long-baseline neutrino oscillation experiments, which 
may be accessible in the present moment. This might 
already give us an opportunity to explore the relevant 
parameter region. 
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III. DIFFUSE NEUTRINO FLUX 


In order to calculate diffuse neutrino flux observed 
at IceCube, we numerically solve the simultaneous partial 
differential equations with respect to differential number 
density (E,^, , z) o f cosmic neutrino Vi, which are given 
in Refs. lIltlMlloi: 


d 


dn^. 

+cncvB 


bn^. + Cvi - cnci/Bhy. ^ vv) 

E 




dCT(r/fei/f^® v,v) 


j,k 




dE^ 


cncuB 


u 


- da{vkv[ 

dEi/j^ 

E^. 


Ci^B 


V^V 


dE,, 


(7) 


d 


dfia- _ 

~w ~m. 

+cncuB 


&hp. + £p. - cncvBnvi ^ cr( 


- CuB 


E 

J.fe ' 


cncvB ^ 
j,k 


pOO 

Ep. 

pCO 




dEi,,rit, 


da{vkv^''^ -)• v^v) 
da{vkvf'''^ D^v) 


dE^, 


( 8 ) 


where i,j,k = {1,2,3} are the indices for neutrino 
mass eigenstates. The time t is related to redshift z 
^ = ~(1 + z)H{z). Following the discussion in 
Ref. [l09j |. we treat cosmic neutrino as an incoherent 
sum of mass eigenstates. The first term on the right- 
hand side is responsible for energy loss of cosmic neu¬ 
trino, owing to redshift, and the energy-loss rate b is 
given with b = E[{z)Eu. The second term represents 
the influx from sources of cosmic neutrino. In this study, 
we assume that all sources provide the same spectrum of 
cosmic neutrino, i.e., Cv^{E^^,z) is simply parametrized 
as L^^{E^^,z) = W{z)Co{E^.) with the cosmic neutrino 
spectrum Co{E^) from each source and the source distri¬ 
bution W{z) with respect to redshift 2 . Here, the source 
distribution function is assumed to be common for all 
the mass eigenstates of cosmic neutrino. We adopt a 
power-law spectrum, which is characterised by the spec¬ 
tral index Si, and the cutoff energy Ecut- 


Co{E,) 


QqE^ exp 


E, 

Ecut 


(9) 


where Qq is the normalization of the flux, which will 
be adjusted so as to fit to the observed flux. This type 
of spectrum typically results from hadronuclear process 
{pp inelastic scattering) in the cosmic-ray reservoir, and 
the values of s^, and Ecut are expected to be determined 
by p rope rties (acceleration rate, i.e., magnetic field and 
size |llOl| ') of the cosmic neutrino source. The flavour 
composition of cosmic neutrino from pp reaction is ex¬ 
pected to be (r'e, iZp, Vr', Eg, Ur) =(1,2,0;1,2,0) at 


each source, which leads to each mass eigenstate produc¬ 
ing approximately with an equal rate, i.e., the normaliza¬ 
tion factor Qq is assumed to be common to all the mass 
eigenstates in our calculations. Although the sources of 
the PeV cosmic neutrinos have not been identified yet, 
we assume the follow ing f unction inspired by the star 
formation rate (SFR) |lll| as a test distribution: 


yv{z) 


(l-fz)3-4 0<Z<1, 

(1 + 1 < .Z < 4. 


( 10 ) 


The third term represents the outflows caused by the 
scattering process with CuB. Here, the cross section cr is 


a{uiu^''^ —>• uu) 


\9'^\^9h 

6tt (s-A/|,) 2-pM|,r|, ’ 

( 11 ) 


where Tz' = PziMz' is the total decay width of 

Z'. The coupling gE = gz'{E^)iaQapUpj is the coupling 
in the mass eigenbasis, where U is the lepton mixing ma¬ 
trix. The number density of cosmic neutrino background 
is given as ncvB = 56(1 -|- z)^ /cm^ for each degree of 
freedom. The constant c appearing in the third, forth, 
and fifth terms is the light speed. The forth and the 
fifth terms provide the influx from the final states of the 
scattering process, the so-called regeneration terms. The 
differential cross sections for cosmic neutrino Ui are cal¬ 
culated to be 


dcrjukuy’'^ Uju) Ig'fcPS; Wu? rn^^E.^ 


dE, 


2tt 


Vk 




, ( 12 ) 


da{ukuf^^ u,u) _ Y.I Wu? m^. {E^^ - 
dE^ 


27r 


^2 


(s-M|,)2 + M|,r2^/ 


(13) 


For cosmic antineutrino Ui, the differential cross section 
for the scattering process Ukv?"^ kiU {ukuf''^ Uiu) 
is the same as Eq. (TT^ (Eq. (fl^ b We numerically solve 
these simultaneous partial differential equations Eqs. dZD 
and ([H]) of cosmic neutrino propa gation, following the al¬ 
gorithms introduced in Ref. |112| | . We confirmed that our 
numerical method correctly reproduces the results given 
in Ref. [4l| . After the simultaneous equations are solved, 
the differential number density of cosmic neutrino at 
the Earth (z = 0) is obtained, and the neutrino flux 
observed at IceCube is calculated as 

cP.AEu^) = En^XE.,.z = d). (14) 


In the next section, we plot the total fluxes 4> = + 

pu-) as functions of the observed energy of cosmic neu¬ 
trino. 
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Lightest neutrino mass (^3) [eV] 


Lightest neutrino mass (m^) [eV] 


Lightest neutrino mass (m^) [eV] 


FIG. 2. Resonance energies as functions of the lightest neutrino mass. The mass Mz' of the leptonic gauge boson is taken to 
be {10, 20,100} MeV, and the redshift z of the scattering point is varied between z = 0 (solid) and z = 2 (dashed). The band 
labelled with ‘Vi” {i = {1, 2, 3}) is the region of the resonance energy corresponding to the scattering with the mass eigenstate 
Vi of cosmic neutrino background. The vertical lines with T „(2 = {0, 2}) indicate the temperatures of CvB at 2 = (0, 2}. All 
the neutrino masses in our reference choices, which are also marked with the vertical lines, are significantly larger than the 
temperatures; The inclusion of the thermal distribution effect of Cz^B momentum does not change our numerical results. 


IV. NUMERICAL RESULTS 

We numerically solve Eqs. 0 and ([ 8 ]) and calculate 
the cosmic neutrino flux in the presence of the — Lr 
interaction. Throughout the numerical study, we do 
not take account of the effect of thermal distribution of 
the Ci^B momentum. This treatment is justified by our 
choices of the parameters; we are mainly interested in 
the parameter region where the lightest neutrino mass is 
much larger than the Cz/B temperature, cf. Fig. [2j We 
have checked numerically that the Cz^B momentum effect 
does not drastically change the conclusions drawn in this 
section. We come back to this point at the end of Sec. 
IV-A. 

In the following calculations, we use the best-fit values 
of the mi xing angles and the mass squared differences 
from Ref. |ll3l| : 

sin^ 6 »i 3 = 0.0234 (0.0240), sin^ 623 = 0.567 (0.573), 
sin^ 6 »i 2 = 0.323, Amj^ = 7.60 x 10“*^ [eV^], 

|Am|J = 2.48 (2.38) x lO’^ [eV^] (15) 

for the normal (inverted) mass hierarchy, and the CP 


violating Dirac phase is set to zero. 


A. Reproduction of the gap with the SFR 

To begin with, we check if the gap can be reproduced 
by the resonant — Lr scatterings. Since neutrino con¬ 
sists of three generations, we generally expect three gaps 
in the spectrum. In Fig. [21 the resonance energies E^es are 
plotted as functions of the lightest neutrino mass wiightest 
with different values of Mz' and with both the mass hi¬ 
erarchies. It can be read off from the plots that in order 
to reproduce a single gap, the masses of neutrinos should 
be quasidegenerate so that the three resonance energies 
are located at the same point. In Fig. |3l we show the 
cosmic neutrino flux with the attenuation effect by the 
Lfj_ — Lr force and compare the results with three different 
values of the spectral index s^. Here we take the normal 
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FIG. 3. The cosmic neutrino fluxes calculated with the 
T/i — Lt gauge interaction are compared with the three- 
year IceCube data [^. The model parameters are taken as 
Mz! = 11 MeV and gz' = 5 x 10“^. The lightest neutrino 
mass is set to be mi = 0.08 eV and the normal mass hierarchy 
is chosen. The SFR is assumed as the redshift distribution of 
the cosmic neutrino sources. The cutoff energy of the original 
flux is placed at Ecut = 10^ GeV. The three different values 
of the spectral index Sv are examined. 


hierarchy with the lightest neutrino mass mi = 0.08 eMl 
and set the model parameters as Mz' = 11 MeV and 
gz' = 5 X 10“"^. For the sources of cosmic neutrinos, 
we assume the SFR, which is given in Eq. CHD, as their 
redshift distribution, and the cutoff energy Ecut, which 
appears in Eq. ([51), is taken as Ecut = 10^ GeV. The 
normalization factor Qq is adjusted so that the magni¬ 
tude of the calculated flux fits the observation. As can 
be seen from the figure, the flux is significantly atten¬ 
uated around 400 TeV — 1 PeV. With a spectrum in¬ 
cluding the gap, one can expect a relatively good fit to 
the observation, although the gap will be shallower than 
the bottom of the calculated spectra once the curves are 
averaged over each energy bin. Since the spectrum calcu¬ 
lated with the inverted hierarchy is essentially the same 
as the normal hierarchy shown at Fig.[3l we do not repeat 
it. 

Let us mention the possibility of simultaneous repro¬ 
duction of the gap and the edge. In view of Refs, [ifll - ll^ . 
we here take lower values of and try to form the edge 
at the upper end of the spectrum by means of the — Lt 
interaction, instead of setting the cutoff energy by hand. 
Note that with an appropriate adjustment of the flux 
normalization, lower values of the spectral index can still 
give a good fit to the current observed spectrum 01 . 
According to Fig. [2l the mass of the lightest neutrino 


This leads to — 0.25 eV, which is slightly higher than 

the 95% C.L . from the combined analysis of cosmological obser¬ 
vations El. However, once the cosmological model is extended 
to include more parameters, the constraint is expected to be re¬ 
laxed. For instance, simu ltan eous inclusion of A^eff mu 

leads to '^mu < 0.28 eV El. 


FIG. 4. The cosmic neutrino flux calculated with Mz' = 9 
MeV and gz' = 4x 10“^. Here the normal hierarchy is chosen 
and the lightest neutrino mass is set to be mi = 6 x 10“® eV. 
The spectral index is taken to be Si, = 2.3 and 2.1. 


should be smaller than 10“^ eV to split the resonance 
energies and distribute them to the positions of the gap 
and the edge. The mass of Z' should be smaller than 
^z' ^ 20 MeV to place the resonance energies at the 
appropriate positions, cf. Eq. (g)). In Fig. 0] ([S]), we set 
the mass of Z' to 9 MeV, the coupling gz' to 4x 10“^, and 
the lightest neutrino mass mi (m 3 ) to 6 x 10“^ eV with 
the normal (inverted) hierarchy of neutrino mass. Here, 
the cutoff energy is taken to be sufficiently high so that 
the numerical results do not depend on the value. The 
gap is successfully reproduced by the scattering with the 
heaviest mass eigenestate of Ci^B. On the other hand, 
the resonant scattering for the edge seems insufficient: 
the flux is attenuated only between 3 and 7 PeV, which 
may be too narrow (and also too shallow) to explain the 
required property of the edge, although it is consistent 
with the current data. 

Lastly, we comment on the effect of the CvB momen¬ 
tum. If the lightest neutrino mass is chosen to be as light 
as the CizB temperature, the CizB momentum effect is 
expected to become appreciable, which would make the 
width of the edge wider. We will study this possibility in 
the near future. 


B. Source distributions 


So far, we have adopted the SFR as the redshift dis¬ 
tribution of cosmic neutrino sources in our calculations. 
However, the source has not been specified yet, and some 
of the ast rophysi cal objects have been discussed as the 
candidate pH - Hq . In Fig. [51 w e exa mine the distribution 
of gamma-ray bursts (GRBs) [ll5j| . 


>Vgrb(2) oc 


(1 + 2 ) 4-8 
(I+ 2)4-4 


0 < 2 < 1 , 

1 < 2 < 4.5 


(16) 

























7 



FIG. 5. The same as Fig. |4] with the inverted mass hierarchy, 
and the monotonic distribution [m, 

Wmono(^)0c(l+2r (17) 

with m = 2, Zmax = 4 and to = 5, Zmax = 1, and com¬ 
pare them to the result calculated with the SFR. Here, 
the spectral index is set to be = 2.3. The neutrino 
mass spectrum and the model parameters are taken to 
be the same as Fig. [3l It is natural to conclude that, 
if the power-law spectrum Eq. (I9|) is assumed, the types 
of source distribution do not make a big impact on the 
shape of the flux. One can expect to find the difference 
caused by the choice of source distribution at the energy 
region slightly below the resonance point, at which the 
effect of regeneration is relevant. 


V. DISCUSSION AND CONCLUSIONS 

We have introduced an anomaly-free leptonic force me¬ 
diated by the gauge boson with a mass of the MeV scale 
in order to simultaneously explain the two phenomena 
with different energy scales in lepton physics: (i) the dis¬ 
agreement between experimental measurement and theo¬ 
retical predictions in muon anomalous magnetic moment, 
and (ii) the characteristic features of the cosmic neutrino 
spectrum reported by the IceCube Collaboration. Assum¬ 
ing that the PeV cosmic neutrinos are produced after the 
pp inelastic scattering process in cosmic-ray reservoirs, 
we have calculated diffuse neutrino flux with the new 
leptonic force. 

We have discussed the relevant constraints, such as 
the lepton trident process and the observation of a solar 
neutrino event at the Borexino, and scanned the model 
parameter space. We have found the choices of parame¬ 
ters, which successfully reproduce the measured value of 
muon anomalous magnetic moment and the gap between 
400 TeV and 1 PeV in the IceCube spectrum. 

Setting the mass of the leptonic gauge boson to be 
around 10 MeV and the lightest neutrino mass to be 
heavier than 1 x 10“^ eV, we could arrange the three res¬ 



FIG. 6. The same as Fig. [3] with different types of source 
distributions. The spectral index is taken as = 2.3. 


onant energies -Eres corresponding to three mass eigen¬ 
states of cosmic neutrino background to the energy 
ranges of the gap = 400 TeV-1 PeV) and the edge 
{El, ~ 3 PeV) simultaneously. However, the resonance at 
the energy corresponding to the edge might be too nar¬ 
row (and too shallow) to explain the sharp upper end of 
the spectrum, which is expected from the observation of 
the IceCube. If one considers the parameter region where 
the lightest neutrino mass is much lighter than 0(10“^) 
eV, momentum distribution of begins to have an 
impact, and one can expect that the inclusion of the ef¬ 
fect would make the resonant region wider. We are now 
preparing for our next analysis in which we examine the 
thermal effect of with smaller neutrino masses. 

We have also examined the different redshift distri¬ 
bution of cosmic neutrino sources and compared the re¬ 
sulting spectra with each other, assuming the power-law 
spectrum for the original flux. One can expect the visible 
difference in the spectrum just below the energy region 
of the gap. 

Before closing this paper, we make some comments on 
Refs, [i^ and [i^l in which the authors also considered 
the Lfi — Lr model and discussed its impact on the cos¬ 
mic neutrino spectrum as well as muon anomalous mag¬ 
netic moment. In Ref. (d^ . the authors assumed coherent 
propagation for cosmic neutrino and calculated the sur¬ 
vival rate of neutrino flavour states for the resonant scat¬ 
tering by considering the evolution of the density matrix. 
In contrast, we assume that parent particles, i.e., cosmic 
rays, are disturbed by the environment and do not decay 
into “free space” at the source of cosmic neutrino, which 
leads us to the equation of incoherent (mass eigenstate) 
propagation. With the cross sections Eqs. (IIIl)-(ini) in 
mass eigenbasis, we calculate the diffuse neutrino flux 
by solving the neutrino propagation equations Eqs. © 
and ([S]). Namely, our treatment of cosmic neutrino and 
approach are different from those adopted in Ref. (d^ . 
Both the treatments of neutrino propagation reproduced 
the gap in the cosmic neutrino spectrum. However, the 
flavour ratios resulting from the two methods are distinct 


























from each other: More electron neutrino can survive with 
the coherent propagation. In addition, we investigate the 
dependence of the spectral index and the source distri¬ 
bution on the flux. In Ref. [i^, the same propagation 
equations are used, and the diffuse neutrino flux is calcu¬ 
lated for several values of neutrino masses and a Z' mass. 
Among their parameter choices, the authors paid special 
attention to the case where the scattering effect becomes 
significant. To illustrate it, they set Mz’ — 1 MeV and 
^lightest ^ 10“^ eV, where the effect of C^B temperature 
made the resonant region wide. On the other hand, we 
choose Mz' — 10 MeV and Slightest ^ 10“^ eV, which 
is tailor made for the reproduction of the gap and the 
edge. With this parameter choice, where the Ci^B tem¬ 
perature effect is not essential at all, we have further in¬ 
vestigated the dependence of the spectral index and the 
source distribution on the flux. In this sense, our study 
is complementary with Ref. (47l |. Furthermore, we have 
discussed the experimental bounds on the leptonic force 
and included the constraints from the loop-induced ly-e 
scattering process [o^ and the BBN [i^, both of which 
disfavour the possibility of the Z' being as light as 0(1) 
MeV. We have confirmed that our numerical method cor¬ 
rectly reproduced the results of Ref. [l^ in the parameter 


region with ruiightest >10 ^ eV. 
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